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Abstract. I argue that the cold dark matter (CDM) model requires that even within a few 
kpc of the center of a galactic halo a significant fraction (greater than a few percent) of the 
surface density is contained in substructures with masses ^ 10^ M©. These structures should be 
light enough to avoid dynamical friction and dense enough to avoid tidal disruption. I then show 
using the results of numerical simulations that this substructure will significantly alter the flux 
ratios of multiply imaged quasars (QSOs) without changing the image positions. The degree to 
which this occurs will depend on the angular size of the QSO and thus the wavelength of the 
observations. 



1 Introduction 

The ACDM model of structure formation has so far been enormously successful in explaining the 
large scale structure of the universe and the anisotropies in the cosmic microwave background (CMB). 
However the model appears to face some difficulties explaining observations on galactic and smaller 
scales. A number of such problems have been discussed in the literature [13, ^, |2|, |l5|, for example]. For 



the present discussion the most relevant of these problems is the observation that CDM simulations 
of the local group of galaxies predict an order of magnitude more dwarf galaxy halos with masses 
greater than ~ 10^ Mq than are observed ^, ^. This could be a sign that there is something 
fundamentally wrong with the CDM model [yTH, %■ Alternatively, the small DM clumps could exist, 
but not contain observable dwarf galaxies. This situation can easily, perhaps inevitably, come about 
through the action of feedback processes in the early universe. For example, photoionization can 
prevent gas from cooling and thus inhibit star formation in small halos |^] . In section |2|, I show that 
the overabundance of DM clumps with respect to visible galactic satellites is likely to extend down to 
smaller masses and larger fractions of the halo mass than have thus far been accessible to numerical 
simulations. 

These nearly pure dark matter structures have largely been considered undetectable (see for a 
review of possible methods). Acting as individual gravitational lenses clumps with velocity dispersions 
of 10 km s~^ (corresponding to a mass of order 10^ Mq) rarely create multiple images and when they 
do the image separations are too small to be resolved (milli-arcseconds) and are not variable on an 
observable time scale. However, as shown in [|lO|, if the CDM model is correct and these substructures 
exist within the lenses responsible for multiply imaged QSOs they will have a dramatic effect on the 
image magnifications. This is the subject of section 0. 



2 The abundance of small mass substructure 

In the CDM model galaxy halos are built through the merger of smaller halos. The larger of these 
subclumps (^ 10^ Mq) sink into the center of the halo through the action of dynamical friction where 
they are destroyed by tidal forces. The orbits of smaller clumps decay much less quickly, but they 
will suffer tidal stripping and collisions. The Nbody simulations of Klypin et al. 0] and Moore et al. 
[^] have shown that substructures of masses ^ 10^ Mq do survive in significant numbers in galactic 
halos. Within the virial radius of a galactic halo roughly 10% to 15% of the mass is contained in such 
subclumps. Here I consider how many smaller mass subclumps will survive. 




Figure 1: The penetration depth of a subclump as a function of the collapse redshift of the subclump, Zcoii, and the 
collapse redshift of the host halo, z^^l° . The scale size of the host halo is Rs (about 20 kpc for the Milky Way). 



Navarro, Frenk &: White ||T^ (NFW) found that CDM clumps in their simulations have universal 
spherically averaged profiles of the form 
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Once the cosmological parameters and the initial spectrum of density fluctuations have been chosen 
the NFW profile is approximately a one parameter family; we can consider r<j to be a function of 5s- 
There is a significant scatter about this relation however. In agreement with the spherical collapse 
model the central density is found to be proportional to the average density of the universe at the 
time when the clump became nonlinear, the collapse time, so 5^ tx (1 + Zcoii)'^ ■ 

Clumps orbiting within a larger halo will have mass stripped off them until their radius reduces 
to approximately the tidal radius, defined by 
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where M{R) is the mass of the host halo interior to R and m{r) is the same for the subclump. The 
superscript ^^halo" refers to the host halo. It was shown in |]lO[ that while the tidal radius is larger 
than the subclump scale length, r^, tidal stripping does not remove a significant amount of mass from 
the subclumps. With this in mind a penetration depth, Rp, can be defined by the relation 
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At galactic radii R <^ Rp mass is transfered between the substructures and a smoother component 
of DM. Since subclumps are not generally on circular orbits Rp should be thought of as the smallest 
pericenter distance a clump can have during its orbits without losing much of its mass. It turns out 
that Rp can be expressed in terms of just 5s/ 5^°"^" = (1 + Zcoiif" /{^ + z^oifY- The penetration depth 
is plotted in figure |l[ It can be seen there that if (1 + Zi.oii)/{'i- + z^oli) ~ 3 or larger the subclump 
can remain intact down to within a few kpc from the center of the host halo. If we are interested in a 
galactic halo that is assembled at z'^^lf ~ 1 then a small mass clump that formed at z^ou ^ 6 is likely 
to survive within it. 

How plentiful should such subclumps be? The merger history of a halo can be calculated using 
the extended Press-Schechter formalism Figure ^ shows what fraction of the mass in a typical 
galactic halo at z = 0.6 is in separate clumps of a given mass range at higher redshifts. The collapse 
of a typical galaxy halo of mass ~ 10^'^ becomes nonlinear at z ~ 1 in the i^matter = 0.3, Vt\ = 0.7 
model used here. Fifty percent of this halo is in clumps of mass m <^ 10^ at z = 6 which must 
have undergone collapse at some earlier time (typically this mass scale collapses at z ~ 7 though in 
the overdense region surrounding the future large halo the mean z^aii could be larger). Even at z = 2, 
15% of the final halo is in free floating clumps of m ^ lO"^ which typically collapsed before z ~ 6. 
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Figure 2: This plot concerns the distribution of progenitors of a halo that is 2 x 10^^ Mq at z — 0.6. The given fraction 
(25%, 50% and 75%) of the halo mass is contained in clumps of masses smaller than the curve at the given redshift. For 
example, 75% of the mass that will end up in the halo is contained in separate clumps with masses < 10^ M0 at z — 10. 
Plotted as a dotted curve is the typical mass of objects undergoing nonlinear collapse. 



It is difficult to imagine how all this very compact substructure could be efficiently destroyed 
within the final halo. Collisions are not effective because they are rare and generally at too high a 
speed to disrupt the clumps. Even when (^) is violated the lose of mass in substructure is a steeper, 
but still gradual function of R. In addition, the observational constraints on dark satellites of the 
Milky Way are rather weak in this mass range |10|. As a result I think it is reasonable to assume 
that at least several percent, perhaps significantly more, of the surface density at a projected radius 
of i? ~ kpc is contained in substructure. A firmer prediction will ultimately require simulations with 
mass resolutions several orders of magnitude smaller than are now available. 



3 Compound gravitational lensing 

Although the gravitational pull on a light beam from a m ^ 10^ Mq DM clump is very weak, if the 
beam also passes through a larger mass concentration the contributions to the magnification from the 
two structures will combine nonlinear ly and enhancing the infiuence of the smaller clump. With the 
larger lens producing multiple images the effects of the small clumps can be seem in the magnification 
ratios. For this to happen the angular size of the source cannot be much larger than that of the 
substructures. QSOs are the obvious choice. 

The large number of subclumps and the small, but not vanishingly small, size of the source make 
numerical simulations of the lensing necessary. Generally an image is infiuenced by several subclumps 
acting together rather than an individual one. As a simple example we consider a singular isothermal 
sphere (SIS) for the primary lens, M{R) oc R. Two images are formed, one farther from the center 
of the lens, image 1, and one closer, image 2. Image 2 is reversed in one dimension with respect to 
image 1 and the original source. The Einstein ring radius which characterizes the image separations 
is re = A'K{ahaio/c)'^DisD\/Ds on the lens plane. Dg; A and Dis are the angular size distances to 
the source, the lens and between the source and lens. The one-dimensional velocity dispersion of the 
host halo is cJhaio- In a flat fimatter = 0.3 model = 3.5 h^^ kpc or 0.3 arcseconds for a lens with 
(Thaio = 150 km s~^ at z = 1 and a source at z = 3. 

In our simulations the subclumps are positioned randomly, but with an average number density 
proportional to the total surface density of the primary lens. They are truncated at the tidal radius 
appropriate for their internal structure and position. This neglects the fact that subclumps may be 
on elliptical orbits which have taken them closer to the center halo in the past. The mass spectrum of 
the clumps is taken to be dN/dm oc in accordance with numerical simulations that extend down 
to 10^ Mq. The lower mass cutoff in the simulations is effectively set by the smallest mass that can 
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Figure 3: The two panels on the left are the magnification distributions for the two images of a QSO at 2: = 3 lensed by 
a singular isothermal sphere at z = 1. On the right is the distribution of magnification ratios expressed in magnitudes, 
Am = 2.5 log(|/ii 1/1/^2 1)- The subclumps have masses of 10'' M0 < m < 10* Mq and constitute 5% of the halo surface 
density. The vertical lines show the values expected for the same lens without substructure. The source is at a projected 
distance of 0.12 Einstein ring radii form the center of the lens. The source is 10 pc in radius. 5,000 realizations are 
shown. 



have a significant effect on tlie lensing, 



m, ~ -4- (4) 
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where Is is the proper size of the source. Because of this cutoff the source size dictates what mass scales 
are being probed by compound lensing. For the same example used before rUc — 10^ M(7){lsh/pc)^^'^ 
at i? = 1 kpc. The broad-line region of a QSO is believed to be a few pc in size. For a more detailed 
description of these numerical simulations see |lO| . 

The simulations show that the centroid positions of the images are not change significantly by the 
substructure unless they happen to lay very near a subclump of mass ~ 10^ Mq or larger which is a 
rare occurrence. Because of this the image positions can still be used to constrain the overall shape 
of the halo. The resulting "smooth model" can then be used to predict the magnification ratios of 
quadruply imaged QSOs. This is a more difficult task than might be assumed because smooth models 
that fit the positions well are usually highly degenerate, predicting a large range of magnification 
ratios. However, progress has been made in this direction which will be reported elsewhere |0, 

Sample magnification and magnification ratio probability distributions are shown in figureM^ These 
are produced by creating random realizations of the subclump positions and masses. Images 1 and 
2 are affected differently by the substructure. This is a result of the one dimensional parity flip 
of image 2. The most likely fi is biased in opposite directions for the two cases which makes the 
magnification ratio biased low; that is closer to equal magnification. Negative values of Am are cases 
where the order of brightness is the opposite of what is expected from the smooth model. The spread 
in Am is quite large compared to typical observational errors, ~ 0.01 mag. In this example Am is 
more than 0.2 mag from the smooth model in 74% of the cases. Note that subclumps contain only 
5% of the surface density in this example. With sufficiently secure smooth models the absence of 
such substructure would be clear in the existing sample of quadruply imaged QSOs if it were not for 
microlensing by ordinary stars. 

There are other approaches to detecting compound lensing that could also be pursued. One is to 
look at the statistics of QSO magnification ratios in general. Substructure will create an overabundance 
of near equally magnified images. Substructures will also deform multiply imaged radio jets on milli- 
arcsecond scales. The images can be comparing for signs of small scale distortions that are inconsistent 



with a smooth lens. These approaches are discussed further in |10|. 



4 Discussion 



The continuum emission region of a QSO in visible wavelengths is very small 1,000 AU) and as 
a result it is susceptible to microlensing by ordinary stars. It has been established that microlensing 
does cause changes in the image brightness of Q2237+0305 as large as a magnitude |17]. This could 



potentially wash out any signal from DM substructures. Fortunately the emission in radio and in 
atomic and molecular lines are believed to come from a larger region. Variability in the QSO could 
also be a complication. The time delays between images range from a fraction of a day to months. 
As long as the source is more then a light-month in size it should not vary on a timescale that would 
interfere with these investigations. These considerations demonstrate the need for magnification ratios 
measured in multiple frequencies. With coming data we may soon have a new and unique test of the 
CDM model. 
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